Electronic Structure of Silicon-Based Nanostructures 
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We have developed an unifying tight-binding Hamiltonian that can account for the electronic 
properties of recently proposed Si-based nanostructures, namely, Si graphene-like sheets and Si 
nanotubes. We considered the sp^'s* and sp^ models up to first- and second-nearest neighbors, 
respectively. Our results show that the Si graphene-like sheets considered here are metals or zero- 
gap semiconductors, and that the corresponding Si nanotubes follow the so-called Hamada's rule 
[Phys. Rev. Lett. 68, 1579 1992]. Comparison to a recent ab initio calculation is made. 
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I. INTRODUCTION 

After the first synthesis of carbon nanotubes (CNT's) 
by lijinia more than a decade ago,^ other types of nan- 
otubes have been predicted and experimentally observed 
such as GaN, BN, and AIN among others. ^"^ However, 
it was not until fairly recently that the most obvious al- 
ternative candidate for creating graphene-like sheets and 
nanotubes was proposed: Si.^ In addition, different struc- 
tures have been proposed, each one with different hy- 
bridizations: sp^ , sp^ — sp^ and sp'^ . So far, most studies 
agree that the sp^ configuration is the least favorable one 
and, in contrast, the sp'^ configuration is one of the most 
favorable ones because of stability reasons.®"^ 

On the experimental side, there are now at least six in- 
dependent reports of the fabrication of silicon nanotubes 
in the laboratory. ^*^^^ On the theoretical side, only a 
handful of papers have explored the electronic properties 
of these nanomaterials; moreover, practically all of them 
correspond to ab initio calculations. ^"^'^^ The above early 
work was reviewed by Perepichka and Rosei.^^ 

In the present work, we apply tight-binding (TB) mod- 
els, so successfully used to study the electronic properties 
of graphene and CNT's, ^^ to the Si nanostructures that 
have sp^ and sp^ hybridization. We will refer to these 
structures as silicene. Si (111), Si hexagonal nanotubes 
(Si h-NT's), and Si gear-like nanotubes (Si g-NT's). Sil- 
icene is a two dimensional sheet with a honeycomb lattice 
of lattice constant a made out of Si atoms which have sp^ 
hybridization. Thus, silicene has the same structure as a 
graphene sheet but it is composed of Si atoms instead of 
C atoms. A Si (111) layer has a lattice structure which is 
the same as the honeycomb lattice for silicene, except 
that one set of atoms (e.g., B) is vertically displaced 
(e.g., down) from the A-plane due to the sp'^ bonding (sec 
Fig. 1). Single- walled Si h-NT's and Si g-NT's arc formed 
by rolling up, respectively, a silicene and a Si (111) sheet. 

Our goal is two-fold. First, we would like to com- 
pare the electronic properties of silicene and Si (111) to 
graphene, and the electronic properties of Si h-NT's and 
Si g-NT's to CNT's. Second, inconsistencies between an 



ab initio calculation and a proposed tt-TB model, both 
by Yang and Ni^ motivated us in developing a coupled 
CT-TT TB model. Results on the effect of the coupling will 
be presented. 

The present paper is organized as follows: in Section II 
we describe the TB theory for the Si nanostructures, in 
Section III we present and discuss our results, and in Sec- 
tion IV we give the conclusions. An appendix has been 
added for further details about the model in question. 



II. THE TIGHT-BINDING MODEL 

We now present a new theory of the band structure of 
a single Si sheet. This theory applies to both the flat Si 
sheet and the Si (111) layer. 

The reason a unifying Hamiltonian is possible is be- 
cause the lattice structure of a Si (111) sheet is simi- 
lar to graphene. Figure 1(a) is a two-dimensional rep- 
resentation of the lattice in question. The A atoms are 
in the ccy-plane and the B atoms are out of it and lo- 
cated at z = — a/(2v6) (the z axis points towards the 
reader). Thus, the word sheet in the present work means 
one atomic plane of A atoms above one atomic plane of 
B atoms, i.e., a sheet is two atomic planes. The ba- 
sis vectors of the lattice are ai = (a/2) (-s/S, — l) and 
3-2 = (a/2) {V^, l) with magnitude a.^^ If we label as I 
the Si-Si bond distance, then a — l/\/2. We will see that 
the choice of this coordinate system facilitates the tran- 
sition from the Si (111) sheet to the silicene one. The 
shaded area corresponds to the two-dimensional unit cell 
of Si (111). Notice that the basis vectors and the unit cell 
of Si (111) are equal to the ones in the honeycomb lattice 
of graphene.^® Due to this similarity with graphene, the 
Brillouin zone of the Si (111) lattice is the same as the 
one of graphene. 

In order to compute the band structure of Si (111), we 
use a first-nearest neighbor (INN) sp^s* and a second- 
nearest neighbor (2NN) sp^ orthogonal tight-binding 
model. For these two models, we derive their respective 
10 X 10 and 8x8 Hamiltonians. These choices were die- 
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FIG. 1. (Color online) Lattice of Si (111), (a) Two-dimensional representation of a Si (111) sheet. The atoms labeled as A are 
all in the xy plane {za = 0) and all the B atoms are located below the plane \z — — a/(2\/6)] . Hence, the sheet is composed by 
two atomic planes: one of A atoms and another of B atoms. The A plane is above the B plane. Notice that the z axis points 
towards the reader. The vectors ai and a2 are the two-dimensional basis vectors and the shaded area is the Si (111) unit cell, 
(b) First-nearest neighbors of Si (111). 



TABLE I. Silicon two-center parameters obtained from Vogl 
et al.^^ and Grosso and Piermarocchi.'^'' Blank spaces corre- 
spond to parameters that do not belong to the model. 

Grosso and Piermarocchi.^" 

-4.0497 

1.0297 

-2.0662 
2.0850 
3.1837 

-0.9488 

0.0000 

0.0000 

0.8900 

-0.3612 



tated by the availability of good Si TB parameters which 
correspond to Vogl et al.^^ [sp^s*) and to Grosso and 
Piermarocchi '^^[sp^). The authors are aware of newer 
TB parametrizations,^^'^^ however, these parametriza- 
tions do not reproduce well the Si bulk bands along the 
KV direction, which is important for the Si nanostruc- 
tures. 

Like silicene and graphene, the wavevector k of Si (111) 
in the Hamiltonian is two dimensional, however, their 
dispersion functions 3j(k), differ due to the a — -k cou- 
pling (see Appendix A). Furthermore, we use the two- 
center approximation (TCA) in order to obtain the an- 
gular dependence in the TB parameters. For the sake 
of completeness, the TCA parameters are reproduced in 
Table L We first specify the position vectors of the INN 



for Si (111): 



Parameter 


Vogl et al.^-* 


Es 


-4.2000 


E^ 


1.7150 


Es- 


6.6850 


{sso)t'' 


-2.0750 


{spa)t'' 


2.4808 


(ppa)t^ 


2.7163 


ippn)^ 


-0.7150 


{s*pcr)AB 


2.3274 


[sso)^^ 




{spa)^^ 




{pp<y)2^ 




ipp^)^^ 





jw = 



fd) 



?(i) 



a a 

a a a 



2^/3'2' 2^6, 
a a a 



(1) 



2\/3' 2' 2^6 



and for the 2NN: 



(0,a,0), 





t(2) 






For simplicity. Fig. 1(b) only shows the position vectors 
of the INN. Notice that the x- and y-components of the 
INN in Si (111) correspond to the ones of graphene. The 
non-zero z-components are the vertical displacements of 
the B atoms of the Si (111) sheet. In other words, the 
choice of the coordinate system facilitates the transition 
from Si (111) and silicene by making the z component of 
the position vectors in Eq. (1) equal to zero. The 2NN 
coincide for both sheets. 

The dispersion i?(k) for the Si (111) sheet is shown in 
Figs. 2(d) and 2(e). For the sp^ model, we can find ana- 



lytic formulae for the dispersion relation at the F point, 



Ep+ (r) — Ep 



3[{ppa)^^ + ippn)^^] 



^[^{PP^K 



AB 



5(pp7r)i 



AB^ 



(3) 



Ep_{r) ^Ep + 3[{ppa)^^ + (ppn)^^] 



--[4(ppcT)r + 5(pp^)n, 

where Ep± is twofold degenerate. 

Next, we explain how this formalism developed for 
Si (111) can be transfered to silicene. 

The TT and a bands in silicene can be obtained from 
the previous Hamiltonian by making the z-component of 
the INN equal to zero and substituting the appropriate 
direct cosines of the INN position vectors. As a result, 
the TT and a bands are decoupled, as in graphene. This 
allow us to consider the bands independently. 

For the tt bands, we recover the well known 2x2 Hamil- 
tonian but now including 2NN interactions. 



H{k) = 



E. 



p < (PP7r)^^.g25(k) 70.912 (k) 

70.912 (k) Ep + ipp7r)^^g25 (k) 



(4) 

where 7o(= Kw")! I) is the transfer integral and the 
(7j(k) functions are given in Eq. (A3). The energy dis- 
persion relation for silicene is found from Eq. (4), 

E{k) =Ep + ipp7T)^^g25{k) ± {pp7T)t''w{k), (5) 



where the dispersion function z«(k) 



l5l2(k)| 



1 -h 4 cos ^%^ cos %2. + 4 cos2 ^ . Notice that if we 

make the 2NN contribution equal to zero, i.e., {ppit)2^ = 
0, we recover graphene's energy dispersion in the INN 
approximation. ^® 

For the a bands, the matrix elements of the Hamilto- 
nian correspond to similar ones given in Eq. (A2) for 
Si (111), but now the gjik) functions correspond to 
the ones given in Eq. (A4). These bands are shown in 
Figs. 2(a) and 2(b). 

By evaluating the Hamiltonian at the F point, we can 
obtain analytical formulae for the dispersion relations, 



^P.±(r)-^p±3(pp7r)f^ 
Es±{T) = E. 



6(ssa)^^±3(ssCT)f^, 



Ep±{T)^Ep + -i[{ppa) 



AA 
2 



AA-\ 
2 



± \ [(PP-)r 



(pPTt) 



where Ep±{^^ is twofold degenerate. 

Finally, the Si h-NT and Si g-NT band structures are 
computed from Si (111) and silicene by imposing bound- 
ary conditions on k along the chiral direction. ^^ The re- 
sults are shown in Fig. 3. 

We now proceed to examine our findings obtained from 
the TB theory. 
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FIG. 2. Band structure of silicene and of Si (111) obtained 
from our TB models compared to the ab-initio results from 
Yang and Ni.^ 



III. RESULTS AND DISCUSSION 

The electronic band structures for the Si sheets and 
the Si-NT's are shown in Figs. 2 and 3. Here, we discuss 
them separately. 



A. Silicene 

The silicene band structure is shown in Figs. 2(a) 
and 2(b) according to the sp^s* and sp^ models, respec- 
tively. Like in graphene, the silicene tt bands are not 
coupled to the a bands due to the planar and orbital 
symmetries.^® When we compare the tt energy bands of 
silicene from the sp^s* model to the corresponding ones 
in graphene, we notice that they have a similar form. 
This occurs because both sheets have the same lattice 
structure. As far as the a bands are concerned, the va- 
lence bands in silicene have been lowered down, therefore, 
the crossings that occur in graphene between the tt and 
a bands do not occur in silicene. Different results are 
obtained when we perform the same comparison using 
the sp^ model. The tt band maintains the form as in 
graphene, nevertheless, the tt* band changes its curva- 
ture as it approaches F from the AT and the MT di- 
rections. This change could occur because of the 2NN 
interactions or because of the signs of the TB parame- 
ters. In order to determine which is the case, one has 
to perform further calculations on the curvature of the 
band, which will be presented somewhere else. In con- 
trast to the sp^s* model, the tt and a valence bands do 
cross. This shows how sensitive the band structure is to 



Present Work 

Si h-NT- (4,4) 



Present Work 

Si g-NT- (4-4) 




FIG. 3. Band structures of Si h- and Si g-NT's according 
to our calculations and to Yang and Ni7 The two-center pa- 
rameters used in our calculations were taken from Grosso and 
Piermarocchi .^^ 



the TB paramctrization of the model. 

It is known that close to the K point, graphene shows 
a linear dispersion, ^■^"^^ 



E{p)=±vM, 



(6) 



where vq and p{= hk) are the Fermi velocity and the car- 
rier momentum, respectively. In graphene, it has been 
measured that vq « 10^ m/s.^'^'^^ The linearity oi E{p) is 
responsible for the electrons to behave as Dirac massless 
fermions.^'^'^^ Figures 2(a) and 2(b) show that silicene 
has a linear dispersion close to the K point. By perform- 
ing a linear fitting, we find that the Fermi velocity vq in 
Eq. (6) is of the order of 10'"' m/s for both models. Com- 
pared to graphene, the electrons move slower in silicene. 
This occurs since the it bonds, which are the responsible 
for conduction in the sheets, are weaker in silicene than in 
graphene (the atomic distance in silicene is greater than 
in graphene). 

Next, let us now compare our sp^s* results to Yang's 
ab-initio ones. The band structure of silicene obtained by 
Yang and Ni^ is shown in Fig. 2(c). When compared one 
notices certain differences. For example, in the neigh- 
borhood of the r point in Fig 2(a), both conduction and 
valence bands have the opposite curvature with respect 
to Fig. 2(c). Possibly, this difference can be attributed 
to the INN approximation in the sp^s* model. There is 
agreement in the silicene n and the lower a bands, how- 
ever, this is not the case for the upper w* band, which 
has opposite curvature. Moreover, there are no crossings 



between the tt and a valence bands. As it was discussed 
previously, this is due to the TB parameters. 



B. Si (111) 

The band structure of Si (111) is shown in Fig. 2(d) 
and 2(e) for the sp^s* and the sp^ models, respectively. 
In Si (111), the sp^ hybridization causes a coupling be- 
tween the TT and the a states. The effect of it is evident 
in the anti-crossing of the originally uncoupled bands in 
graphene. ^^ Note, however, that the w band is still doubly 
degenerate at the K point. In graphene, the degeneracy 
at this point is required by the hexagonal symmetry. ^^ 
When we lower the B atoms in silicene to create the 
Si (111) sheet, as it is shown in Fig. 4, this symmetry 
is not removed hence, the tt bands are still degenerate at 
K in Si (111). 

Notice that the sp^ and the sp^s* models lead to dif- 
ferent bands, particularly, close to the F point. As the 
bands approach this point, the curvature of some of them 
changes. Consider, for instance, the second-lowest va- 
lence band in Figs. 2(d) and 2(e). Close to F, this band 
has negative curvature in the former figure, while it has 
positive curvature in the latter one. These differences 
might be due to 2NN interactions. We computed the 
band structure of Si (111) using two parameter sets for 
a sp^ TB model by Chadi and Cohen;^^ one of them for 
INN only, and the other one for INN plus one 2NN inter- 
action. These results are not presented here, however, we 
mention that the band without the 2NN interaction has 
negative curvature, like Fig. 2(d), and the band with only 
one 2NN interaction has positive curvature, like Fig. 2(e). 
Another possible reason for the curvature differences is 
the sign of the TB parameters. Further calculations on 
the curvature of the bands at the F point are needed in 
order to determine if it is due to the 2NN interactions 
or the TB parameters. These calculations will be shown 
somewhere else. 

As far as the eigenstates are concerned, they are dif- 
ferent in graphene and in Si (111). The it bands in both 
structures are a good example of this. Whereas the eigen- 
states of the the tt bands in graphene are Pz orbitals, the 
eigenstates of the tt bands in Si (111) correspond to a lin- 
ear combination of s,Px,Py and Pz orbitals. This occurs 
because in Si (111), the tt and a bands are coupled. 

Notice in Figs. 2(d) and 2(e) that close to the K point, 
the dispersion relation is linear, which indicates the pres- 
ence of Dirac massless fermions. By performing a linear 
fitting, we find that the Fermi velocity vq is of the order 
of 10^ m/s. Compared to graphene, the Dirac fermions 
move slower in the Si (111) sheet. In order to understand 
the difference in velocities, let us look at the hybridiza- 
tions of both sheets, as it is shown in Fig 4. 

On the one hand, each atom shows sp^ hybridization 
in graphene. In this hybridization, each atom has a lobe 
that is perpendicular to the sheet plane (the sheet cor- 
responds to the xy-plane and the lobe is oriented along 



t^t r^i 



(a) 



(b) 



FIG. 4. (Color online) The sp^ and sp^ hybridizations in (a) 
silicene and in (b) Si (111). 



the positive z-axis as it is shown in Fig. 4(a)). All lobes 
are oriented along the same direction, i.e., the z-axis, 
and, therefore, they form tt bonds with their INN. The 
■K bonds are responsible for the conducting character of 
the sheet. 

On the other hand. Si (111) shows sp^ hybridization. 
Figure 4(b) shows that each atom has a lobe that is per- 
pendicular to the sheet, however, neighboring atoms have 
their lobes pointing in opposite directions, i.e., the pos- 
itive and negative z-axes. Due to the alternating orien- 
tation, a lobe that points along the, say, positive z-axis, 
will not form a tt bonding with its INN, but with its 
2NN. 

Since the tt bonds in graphene occur between INN, the 
coupling is stronger than in Si (111), where the bonds oc- 
cur between 2NN. An electron finds it "easier" to tunnel 
from one atom to the another one when the coupling is 
stronger. This explains why electrons in graphene move 
faster than they do in Si (HI). 

The band structure of the Si (111) sheet computed 
by Yang and Ni.^ is shown in Fig. 2(f). When com- 
pared to our sp^ findings, we find good agreement be- 
tween them, especially, along the KTM directions. Close 
to r, the curvature problem is overcame by including 
some TCA interactions from 2NN in the Hamiltonian, 
namely, (ppu)^'^ and (pp7r)^^. Most of the differences 
occur along the MK direction, for instance, our calcu- 
lations show that the second-lowest valence band should 
have the opposite curvature when it is compared to Yang 
and Ni'' The positive curvature of this band comes from 
the strong repulsion induced by the p-like band located 
above it. 



C. Si h-NT 

Nanotubes are fully characterized by their chiral vector 
C/j = n&i -\- ma2, where n and m are integer numbers, 
and ai and 3.2 are given above. ^* The band structures 
of Si h-NT's with chiral vectors (4,4), (6, 6), (8, 0) and 
(12,0) are shown in the left-most column of Fig. 3. The 
dashed line corresponds to the Fermi level. These bands 
were obtained by substituting the quantized nanotube 



wave vector k in Eq. (5). In this figure, the dispersion 
g25(k), introduced by the 2NN interactions, its hardly no- 
ticeable, therefore, the band structure of Si h-NT's looks 
very similar to the one of CNT's in the INN approxi- 
mation. It is due to this similarity, that we neglect the 
dispersion (?25(k) in the band structure of Si h-NT's and 
discuss it as if it were a INN approximation only. In 
the INN case, the only difference between the CNT and 
the Si h-NT band structures is a scaling factor, which 
corresponds to the transfer integral 70 (70 = —3.033 eV 
for C and 70 = -0.949 eV for Si). Whether or not a 
Si h-NT would be conductor or semiconductor, does not 
depend on 70 but on its symmetries, thus. Si h-NT will 
have similar electronic properties as CNT's, i.e., ii n — m 
is a multiple of 3, the nanotube is a metal, otherwise, 
it is a semiconductor. We will refer to this property as 
Hamada's rulc.^^ From Hamada's rule, we conclude that 
all armchair Si h-NT's are conductors (n ~ m) and zig 
zag tubes are conductors if n is a multiple of 3 {m = 0). 
The band structures for Si h-NT's shown in Fig. 3 con- 
firm this rule. 

We proceed to compare the energy band gaps between 
CNT's and Si h-NT's. According to Dresselhaus et al.^o 
the band gap is given by 



Eg = 27oMini(; fk^ + ^{N±l)Ki 



(7) 



where Ki and K2 are the reciprocal vectors of graphene 
and silicene, 70 is the transfer integral, N is the number of 
hexagons in the nanotube unit cell, k is the wave vector, 
and Min w is the minimum of the graphene and silicene 
dispersion relation w{= |gi2(k)|) with respect to k. In 
particular, for zig zag tubes, the minimum occurs at fc = 
0,^® according to Eq. (7), this yields to. 



Eg = 27oWl + 4cos 



T^qo 



4 cos^ 



(8) 



where qo = (1/3) (TV ± 1) and N = 2n. In Fig. 5, we show 
our results for the band gap of both CNT's and h-SiNT's 
as a function of the diameter. For a given diameter, the 
band gap Eg of Si h-NT's is smaller than the band gap of 
CNT's. In order to understand this result, we should look 
to the approximate expression for Eg derived by Saito et 
al.is and Ando,^* 



En 



\J\_ 
d 



(9) 



where 7 = \/3a7o/2 and d = |C;j|/7r (the nanotube diam- 
eter). For a given diameter. Eg depends only on the pa- 
rameter 7, which depends on the transfer integral (recall 
that the transfer integral corresponds to the interaction 
between neighboring pz orbitals) . Since this interaction 
is smaller for Si h-NT's than for CNT's, the band gap is 
greater for the latter. 

Notice that the 1/d dependence in the nanotube band 
gap is expected. If we let d — >■ cx), then Eg — )■ 0. The limit 
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FIG. 5. (Color online) Band gap of CNT's (circles) and Si 
h-NT's (squares) as a function of their diameters. 
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FIG. 6. (Color online) Effective masses of CNT's (circles) and 
Si h-NT's (squared) as a function of their diameters. 



d — > cxj corresponds to making a nanotube thicker and 
thicker, and, therefore, more similar to graphene, which 
has a zero band gap. 

We proceed to compute and to compare the effective 
masses of Si h-NT's and CNT's. By differentiating the 
band structure of zig zag nanotubes given in Eq. (5), and 
neglecting the 2NN interactions, we obtain an analytical 
formula for their effective masses at the K point. 



dPE/dk^ 



fc=0 



± 



4 cos ■ 



4cos2^:^ 



" 37oa2 



Tqa 



(10) 



Equation (10) is plotted in Fig. 6. First, notice the 
general trend for Si h-NT's and CNT's that the effec- 
tive masses decrease as their diameter increases. This 
occurs because electrons in a Si h-NT (CNT) with 
large diameter, behave more similar to electrons in sil- 
icene (graphene). In other words, if we increase the 
nanotube diameter, we will see that the electrons start 
"losing" their masses, until they become Dirac masslcss 
fermions. Second, notice that for a given diameter, elec- 



trons in Si h-NT's are heavier, in magnitude, than elec- 
trons in CNT's. The top valence band and lowest con- 
duction band are flatter than the corresponding ones in 
CNT's, which explains their mass difference. 

We point out that the effective masses can be obtained 
from k • p theory as well. Ando^'* derived the following 
dispersion relation for the nanotube band gap. 



£;„±(fc)= ±7^^(71)2 +fc'2, 



(11) 



with 



K= — (^n--j,n = 0,±l,±2,... 

where k' is the wave vector measured from the K point, 
L is the magnitude of the chiral vector and \v\ = 1 for 
semiconducting nanotubes. Hence, the effective mass of 
the zig zag tubes is 



(FEQ±/dk'' 



± 



4^2 



3^/37oc 



(12) 



This equation confirms what was stated above about the 
nanotube effective masses using TB theory, that is, m* 
decreases as the d increases since it explicitly depends 
upon d~^. 



D. Si g-NT's 

The second column of Fig. 3 shows the band struc- 
tures of Si g-NT's with chiral vectors (4,4), (6,6), (8,0) 
and (12,0) obtained from the sp^ model. Notice that for 
all graphs but the ones in the second column, the Fermi 
energy is at zero. From this figure, one notices that there 
is a proliferation of the number of bands in the Si g-NT's 
when compared to CNT's. In CNT's there is a total of 
2N bands. ""^^ The factor of 2 comes from the 2x2 Hamil- 
tonian, and N bands come from the quantization of the 
wave vector k along the Ki direction. In Si g-NT's, there 
are 8A^ bands; eight from the 8x8 Hamiltonian, and the 
other N from the same quantization in k. Hence, Si g- 
NT's have four times the number of bands that CNT's 
have. Moreover, there are 4 electrons per Si atom and a 
total of 2N atoms in the g-NT unit cell, hence, there are 
4 X 27V = 8A^ electrons per unit cell. In particular, arm- 
chair and zig zag g-NT's, have 16n bands and electrons 
per unit cell {N = 2n). Table II summarizes some of 
the general characteristics of the h-NT and g-NT energy 
bands. 

After wrapping the Si (111) sheet into a zig zag g- 
NT, we find that the (8,0) tube is a semiconductor with 
a gap of 0.34 eV, and the (12,0) is either a metal or 
a semiconductor of gap zero. The armchair (4,4) and 
(6,6) tubes could be metals or zero gap semiconductors, 
as well. In order to completely characterize the electronic 
behavior of the last three nanotubes, we need their den- 
sity of states (DOS), which is not done here. However, 
due to the similarities with CNT's, we would expect them 



TABLE II. General characteristics of the h-NT and g-NT en- 
ergy bands. Here, e~ and u.c. refer to electrons and the 
nanotube unit cell, respectively. 



h-NT (sp2 



-NT [sp^ 



e~ per atom 


1 


atoms per u.c. 


27V 


Number of bands 


2N 


e~ per u.c. 


2N 



4 
2N 
8N 
8N 



to have a finite DOS at the Fermi level and thus, to be 
metals. 

In general, Si g-NT's obey Hamada's rule like CNT's 
do. In order to understand this, we apply a similar 
reasoning used for CNT's. ^® We first started with the 
Si (111) band structure (c.f. Fig. 2(e)), which has zero 
gap at the K point of its Brillouin zone. Second, we 
sliced the Si (111) bands by quantizing the wavevector k 
along the Ki direction. Third, we used the zone-folding 
approximation to plot the Si g-NT bands. If any of the 
"slices" cuts the K point then, the g-NT is a conductor, 
if it does not, then it is a semiconductor. In the language 
of group theory, the degeneracies caused by the vertical 
and the horizontal mirror plane symmetries in CNT's are 
not removed in the gear-like tube. 

Next, let us compare our Si g-NT energy bands to the 
bands by Yang and Ni.^ From the previous comparisons 
of the Si sheets, we find better agreement between our 
calculations and the ab-initio ones by Yang and Ni'' when 
we use our 2NN sp^ TB model. Hence, we only compare 
our energy bands obtained from this theory to the ab- 
initio ones. 

The Si g-NT energy bands by Yang and Ni^ are repro- 
duced in the left-most column of the Fig 3. We find that 
all our band structures differ from their work. Whereas 
all our armchair tubes are metals or zero-gap semicon- 
ductors, their corresponding tubes are semiconductors. 
Our zig-zag tubes follow Hamada's rule, while theirs do 
not. 

The authors argue that in zig zag Si g-NT's with small 
diameter like (8,0), gaps do not occur because of the a* — 
■K* coupling. We think that this hypothesis is doubtful, 
since our theory does include this coupling and it does 
not close the gap of the small semiconducting tubes like 
(8,0). Furthermore, band gaps tend to open and not to 
close when couplings are added. A possible explanation 
for the discrepancy is that the ab initio calculation suffers 
from the density-functional-theory band-gap problem. 

These significant differences between the TB and the 
ab initio bands are quite surprising since both approaches 
coincide on the band structure of Si (111). We point 
out, though, that in their report, Yang and Ni^ do not 
explain the gap openings when the Si (111) sheet is rolled 
up to form Si g-NT's.^ One could think that curvature 
effects in Yang's results might open a gap in Si g-NT's, 
nonetheless, such gap should be, at most, of the order 



of mcV. Consider a CNT and a Si g-NT with the same 
chirality, for instance, C^ = (12,0). It is known that the 
gap opening /S.E due to curvature in CNT's, is inversely 
proportional to it and it is close to 10 meV.^^''^"'^ The 
transfer integral I7I, can be used as an energy scale for 
the nanotubes, thus, /S.E ^ \^\/(P. The ratio between 
the openings in CNT's and Si g-NT's is then given by 



AS; 



Si 



l7si|/rfi 



ASr 



l7c|/d^ 



^/^Ei 



Si 



iTsil 
l7c| 



rfsi 



l7Si| 
l7c| 



l^Er 



dc 

dsi 



(13) 



Consider a C- and a Si-NT with diameters dc = 
30 A and dsi = 46 A, respectively.^'^® Substituting |7c| = 
3.033 eV, |7si| = 0.949 eV and AEc = 10 meV^'^*^ in 
Eq. (13), we find that 

AEsi ~ ImeV. 

The band gap in the (12, 0) Si g-NT of Yang and Ni'' (see 
band structure in the lower right corner of Fig. 3) is close 
to 200 meV so, curvature cannot be responsible for the 
whole opening of the gap in Yang's nanotubes. 

Our TB scheme, on the other hand, provides an ex- 
planation for the gap behavior in all g-NT's openings in 
terms of Hamada's rule. For example, our (12,0) Si g- 
NT has, according to Hamada's rule, a zero gap since 
this corresponds to a zig zag tube with n multiple of 3. 



IV. CONCLUSIONS 

The electronic properties of silicene. Si (111), Si h- 
NT's and Si g-NT's were studied via a TB approach. We 
derived sp^s* and sp^ Hamiltonians up to INN and 2NN, 
respectively. 

We compared the band structures of Si (111) and of sil- 
icene to the one of graphene. Since all of these materials 
have in-plane symmetry, the n bands are two-fold degen- 
erate at K. We expect Si (111) and silicene to be either 
semiconductors of band gap zero or metals. Electrons in 
the neighborhood of the K point should behave as Dirac 
massless fermions due to the presence of the Dirac cone 
in both structures. However, the Fermi velocities vq hi 
Si (111) (10* m/s) and in silicene (10^ m/s) are smaller 
than the one in graphene (10^ m/s). Electron tunneling 
between Si atoms is less favorable in Si (111) and silicene 
than in graphene because the tt interaction is weaker in 
the the first two cases. 

Silicon h-NT's and Si g-NT's were compared to CNT's, 
as well. The band structure of Si h-NT's and of CNT's 
are similar two each other. Even though we performed 
calculations including INN and 2NN, we found that the 
effect of the latter ones on the bands is negligible. This 
allowed us to make further approximations when calcu- 
lating the band gap and effective masses of Si h-NT's. 
In the case of zig zag semiconductor Si h-NT's, the gap 



is inversely proportional to the tube diameter, as in zig 
zag CNT's, nonetheless, for a given diameter. Si h-NT's 
will have a smaller gap. The magnitude of the effective 
masses is also inversely proportional to their diameter, 
however, for a given diameter Si h-NT's have greater 
mass, which makes CNT's more suitable for transport 
properties. As far as Si g-NT's are concerned, we found 
that they follow Hamada's rule as CNT's do, even though 
they show different hybridizations. 

Our calculations for all the Si-based materials consid- 
ered here were compared to the ab initio calculations per- 
formed by Yang and Ni.^ When comparing silicone and 
Si (111), we found that the 2NN sp^ model is in better 
agreement with Yang's band structures than the INN 
sp^s* model. For this reason, we chose the sp^ model 
in order to reproduce the energy dispersions of Si nan- 
otubes. 

For Si h-NT's, our band structures agree with the ones 
obtained by Yang and Ni:^ they all follow Hamada's rule. 
In contrast, they disagree for all Si g-NT's. Whereas our 
calculations show that these nanotubes also follow this 
rule, Yang's calculations do not. We emphasize that the 
cr-TT coupling does not close the gap for Si g-NT's with 
small diameters, contrary to the hypothesis of Yang and 

I 



Ni.^ Our calculations also show how critical it is to obtain 
accurate TB parameters for applications to Si sheets and 
nanotubes. 
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Appendix A: Tight-binding Hamiltonian 

In this appendix, we explicitly give the Hamiltonian of 
our model. The sp^s* Hamiltonian has the form 



H{k) 



hAA hAB 
hsA hBB 



(Al) 



where the hij 's are 5x5 sub-matrices which arc given as 
follows: 



hAA 



hAB — 



{s/s)aa {s/x)aa {s/y)AA 
{x/x)aa {x/y)AA 

{y/y)AA 

t 







{z/z)aa 



{s/s*)aa 

{x/s*)aa 

{y/s*)AA 



{s*Is*)aa. 



T-BB 



'AB 
AB 



{s/s)ab {s/x)ab {s/y)AB (s/z)ab (s/s*) 

{x/s)ab {x/x)ab {x/y)AB {x/z)ab {x/s*)ab 

{y/s)AB iy/x)AB {y/y)AB (j//z)ab {y/s*)AB 

{z/s)ab {z/x)ab {z/y)AB {z/z)ab {z/s*)ab 

{s*/s)ab {s*/x)aB {s*/y)AB {s*/z)ac {s*/s*)aB 



(A2) 



(sA«)aa -{s/^Taa ~is/y)AA 

ix/x)*Aji {x/y)*AA 

{y/y)AA 

t {^I^Yaa 



{s/s")aa 

i^/s*)AA 

{yh*)AA 

{^^/s*)aa 
(s*/s*)aa 



and hBA = [hAB] ■ In these equations, hAA and hAB and A — B lattice points, respectively. In the case of the 
correspond to the interaction between atoms &t A — A sp'^ Hamiltonian, each sub-matrix is 4 x 4 instead. The 

matrix elements are given as follows: 



{s/s)aa = {s/sYbb =Es + (sscr)^^(7i3(k), 

{s/x)aa = -{s/x)*BB = (spcr)^^.9i4(k), 

{s/y)AA = -{s/vYbb = (spCT)^'^.9i5(k), 

{s/z)aa = -{s/vYbb = 0, 

(a;/x)AA = (a^/a;)BB = -Bp + (ppa)^^.gi6(k) + (pp^)^^gi7(k), 

(a;/y)A^ = (x/y)^s - [{pp<r)2^ - {PPT^)2^] 5i8(k), 

{xIz)aa = {x/z)*BB = 0, 

{y/y)AA = {y/y)*BB =Ep + (ppa)^^.gi9k + (pp^)^^g2o(k), 

(y/z)AA = (y/zYBB = 0, 

(z/z)aa = (^/2)bb ^E.p + (pp7r)^^g25(k), 

(x/a:)AB - (ppa)f ^.g3(k) + {ppir)^'' gi{k) , 

{y/y)AB = (ppa)^^36(k) + (pp7r)f ^<?4(k), 

(z/z)^B = (ppa)^^.gn(k) + (pp^)^^.gi2(k), 

(s/x)ab = -(a;/s)AB = {spa)f^gi{k), 

{s/y)AB = -(y/s)AB = (spcr)f^fl'2(k), 

(s/z)Ai3 = -(2/s)Ai3 = (spcr)f^g8(k), 

(x/y)^B = (y/x)Ai3 = [bpc7)i^^ - (ppTr)^^] g5(k), 

{x/z)ab = (^/2;)as = [(ppfT)f^ - {ppTT)i^] 59(k), 

(y/z)^B = {z/y)AB - [(WT)f^ - (pp^)i^^] 5io(k), 
{s*/s*)aa - (s7s*)bb = E,' + (s*sV)^V3(k), 
(sVa;)^^ - -{s*/xYbb = (s*pa)^^5i4(k), 
(s7y)AA = -{s*/yYBB = (.s>a)^^gi5(k), 
(sVz)^^ - ^{s*/yYBB - 0, 
(s/s*)aa = {s/s*Ybb = (ssV)^^.gi3(k), 
(s7x)^B = -{xIs*)ab = (s*pa)f%(k), 

(s72/)aB = -{yls*)AB = (sV)f^ff2(k), 

(s77^s = -{zIs*)ab = (s7a)rg8(k). 



I 

where the 5y(k) functions are given, respectively, for 
Si (111) and silicone in Eqs. (A3) and (A4). 
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54 (k) = i (e*k.J<^' + Te*"^-^^'' + Te'i'-^^ 

.a(k) = ^(e-^^%e-'^^^'), (A3) 

5s(k) = -i(e-'^i%e-^^%e-'^^^' 



59(k) = — l-(^e-"^ __e-". --e 

^-(p'k-^i +P»k-^2 4-P»k')3 



5ii(k) = - |e-"i +e"""^ + e" 

5'i2(k) = - (e**^"! +e''^''2 + e'" "^ = - e vt? + 2e ^v/a cos -^ 

5i3(k)=e*''-'^i +e'*'-^^ + e'*''^^ + e*''-^* + e'*''^-'- + e'^'^'^ , 
/, \ \/3 / jk-^<" jk-^"' ik-^<" jk-^<^' 



gi5(k)=p. "i +e'''"2 _|__e*k03 ^ _ g^k o^ _|^ _ g»k o^ ^ik o^ 

/, N 3 / jk-J<" lk-^<" tk-^<" ik-^'" 

gi6(k) = - (e*'^"3 +6**^"* +6'""= +e*'^''6 

gij(k) =6'^°^ +6'" "2 +_e*k03 +_e*'^<'4 _^_e»k05 ^ _ g^k o^ ^ 

/, N ^ /" ik-5<" ik-<5'" ik-^'" lk-^<^> 

5lg(k) = — (-e'*'''3 _e»k04 ^^g'kOj ^gikOe 

5i9(k) =e"'''i +e"'°2 + - e '' + - e**^ "-i + - e'" ''f' + - e'" "« 
3/,:k.5f' .k.^f' ,k.^f' ,k.^r 



(2) 



52o(k) - - (^e"'-"^' +e*k.o, + g^k-o, ^g 

52i(k) = g22(k) = g23(k) = 0, 

.g24(k) - e*-'^i" + e'*-^^" + e^'^-^'^' + e^'^^' + e'^^^'^^" + e^"^'^' 



where the position vectors 5, are given in Eqs. (1) and 
(2). 
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50 (k) 

51 (k) 

52 (k) 

53 (k) 

54 (k) 

55 (k) 

56 (k) 

57(k) 

58 (k) 

512 (k) 

513 (k) 

514 (k) 

515 (k) 

516 (k) 

517 (k) 

518 (k) 

519 (k) 

520 (k) 

521 (k) 
525 (k) 



2' 2 



V3 / ,k.5^^' ,k.^<^' 



I 



Jfe-^^Ve-^"' 



^ .^-S^' ^ i ,^-S^^' ^ i .^^s 



(1) 
4 



= 59(k) = 5io(k) = 5ii(k) == 0, 

^ ^.k.^i^> ^ ^,k.<5<^> + ^,k.<5<^' ^ ^.^ ^ 2e-'IS cos ^, 



V3 



'gik-^f ' _ g^k-^f ^ p,k.^<' 



<^' g.k.^f _ ^^k-sr 



■1 A(^) 



2 ^2 ^2 2 



(2) ., r(2) 



4 V 

■1 A(^) 1 /S(^) 1 1 A<2> 1 -1 A<2) 1 -1 A(2' 1 -1 AC 



V3 / g.k.^f _ g.k.Jf 
4 V 



■I A(2' -1 A<2) 



.k.<5f 



„Jk 



^f > ^ 1 „^k.^'^> ^ 1 .k-df' ^ 1 „»k.(5l,=' ^ 1 ,k.^<^ 



3 / -, A'^' -1 A<^' 1 A'^' 1 A'^' 
-le 3 + e "* + e "= +e e 

4 V 

522 (k) = 523 (k) = 524 (k) =: 0, 

jk.^(2) j,,-^'^' ik-^'" 
:g«kOi _|_g»k02 +e'''"3 _)_e 



k.^r , ,k.<5f > , ,k.<5^^' , rk-Sf , ^,k.^f' , .k.,5r> 



(A4) 



r 
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